Nonequilibrium Electron Distribution in Presence of Kondo Impurities 
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We study the energy relaxation of quasiparticles in voltage biased mesoscopic wires in presence of 
magnetic impurities. The renormalization of the exchange interaction of Kondo impurities coupled 
to conduction electrons is extended to the case of a nonequilibrium electron distribution, which is 
determined self-consistently from a Boltzmann equation with a collision term due to Kondo impurity 
mediated electron-electron scattering. The approach leads to predictions in quantitative agreement 
with recent experiments by Pothier et al. [Phys. Rev. Lett. 79, 3490 (1997)]. 
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In 1997 the mesoscopic physics community was puz- 
zled by two experimental findings. On the one hand Mo- 
hanty, Jariwala, and Webb |ffl] have extracted the electron 
dephasing time from weak localization measurements of 
one-dimensional gold wires and found larger dephasing 
rates than predicted by the standard theory of Altshuler 
and Aronov (AA) 0. In the same year, Pothier et al. 
H published data on the nonequilibrium electron dis- 
tribution in mesoscopic copper wires in presence of an 
applied voltage. The electron-electron scattering rate 
gained from these data was shown to exceed predictions 
based on the AA theory. It was intuitively clear that 
these two observations are very likely due to a common 
origin, and in the last three years a large body of work has 
proposed all kinds of mechanisms ranging from interac- 
tion with two-level systems , heating by radiation || , 
1/f noise [0 to intrinsic dephasing at zero temperature 
0. None of these predictions could give a quantitative 
description of the experiments or they were ruled out 
by subsequent experimental studies. Finally, a scenario 
based on the two channel Kondo effect of symmetrical 
two level systems has recently been proposed by Kroha 
H . For a rather complete discussion and analysis of avail- 
able experimental data on both the weak localization 
and nonequilibrium electron distribution measurements 
in various metals we refer to the thesis by Pierre || . 

Already Mohanty et al. jl| have demonstrated that iron 
impurities added to the gold wires can lead to similar ef- 
fects than those observed in nominally pure wires. In 
fact, equating the pair breaking rate in superconductors 



containing magnetic impurities |10| with the magnetic 
contribution to the dephasing rate in normal metals [[ll| , 
one obtains a satisfactory explanation of the weak lo- 
calization data near and above the Kondo temperature, 
where the theory is valid 0. On the other hand, the 
effect of Kondo impurities on electron-electron scatter- 
ing has been addressed only very recently. Kaminski and 
Glazman (KG) |r3] have determined the two-particle t- 
matrix mediated by magnetic impurities for an equilib- 
rium system. The renormalization of the exchange cou- 
pling was studied with poor man's scaling using the ap- 



plied voltage V as a low-energy cutoff. While the the- 
ory can explain the order of magnitude of the observed 
electron-electron scattering rate, it does not reproduce 
the correct voltage dependence but leads to more than 
an order of magnitude deviations from the data over the 
range of parameters investigated experimentally []l3|| . In 
this work we demonstrate that it is essential to go be- 
yond poor man's scaling since deviations from AA pre- 
dictions are pronounced only below the Kondo tempera- 
ture. Most importantly, the renormalization flow has to 
be determined in presence of a nonequilibrium electron 
distribution calculated self-consistently from the Boltz- 
mann equation. The collision term then depends on the 
distribution function not only through the occupation 
probabilities of in- and outgoing electrons but also via 
the renormalized interaction kernel. 

We start by describing briefly the experimental situa- 
tion: A mesoscopic wire of length L with diffusion con- 
stant D is attached to leads biased by a voltage V. (For 
details cf. Refs. ||] and 13 1). The nonequilibrium elec- 
tron distribution function f(e,x) in presence of a steady 
state current can be determined from the Boltzmann 
equation for a diffusive mesoscopic wire [pL4j] 



1 d 2 f(e,x) 



TD 



dx 2 



(1) 



with the boundary conditions /(e, 0) = fp{e — eV/2) and 
f(e, 1) = fp(e + eV/2) imposed by the leads. Here x is 
the position within the wire measured in units of L, and 
to = L 2 1 D. In the simple case of vanishing interaction 
IcoW = we obtain a double step function 

/o(e, x) = (1 - x)f F (e - eV/2) + xf F (e + eV/2). (2) 

The interaction smears these steps, and in the limit 
of strong electron-electron interaction leads to a Fermi 
function with an effective temperature of the electrons 
]l5[ . The distribution function f{e,x) is determined in 
the experiment at various locations by tunneling spec- 
troscopy. 

Assuming a local interaction, the collision integral 
reads 
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/coll (x, e,{/})= /*" / <&K(ui,e,<?) (3) 

x{/(e)/(e')[l-/(e-o;)][l-/(e'+a;)] 
-[l-/(e)][l-/(e')]/(e-o;)/(e' + a;)} 

with an interaction kernel if (w, e, e'). Not to overload the 
notation we suppress here and in the sequel the spatial 
dependence of the kernel and the distribution function. 
To determine the kernel K(uj,e,e'), we follow KG and 
start from the s-d exchange Hamiltonian 



whereby 



H = H + Hr 



H =22 e k Cl a C k 



(4) 



(5) 



describes free quasiparticles with one-particle energies e k 
and creation (annihilation) operators C k (C ka ) of states 
ka. We assume that the density of impurities is small 
enough to treat the interaction with each impurity inde- 
pendently. Then, for a single impurity 



H T = J J2 S-s**'ClC k < 



k'a' 



(6) 



kk' ' aa 1 



where S is the impurity spin operator and s the vector 
of Pauli matrices. Further J is the exchange interaction. 

Let us first address the electron-electron interaction 
mediated by Kondo impurities in an equilibrium metal. 
We rewrite the most singular parts of the interaction in 
terms of single-particle t-matrices and include renormal- 
ization effects for lower temperatures by an approach due 
to Zittartz and Muller-Hartmann Jl6|Jr^| based on the 
Nagaoka equations This theory, though not able 

to describe the zero temperature limit correctly, leads to 
meaningful results at higher temperatures down to tem- 
peratures well below the Kondo temperature jnj . 

For an effective two-particle interaction, the kernel 
K(u>, e, e') is essentially given by the modulus squared 
of the on-shell two-particle t-matrix. Further, one has 
to sum over all final electron spins 07, a'j and the initial 
spin a' of the second electron and average over the initial 
spin a of the first electron and the impurity spin S. For 
impurities with density C lmp the kernel then takes the 
form 



K(uj,e,e') = C i 



E 



(^ll^fe^feV'^fc/oy,*;' a > \ 2 \S) 



aa fa' a', S 



(7) 



where p is the electronic density of states at the Fermi 
level. The two-particle t-matrix is defined by 

Tka,k'a'^k f a f ,k' f a' f = (kf(Tf,k'fCr'f\T\kcr,k'cr') (8) 



and does not depend on the directions of the incom- 
ing and outgoing electrons. Hence, the wave numbers 
k, k 1 and kf, k' f are characterized by the energies e k = e, 
6k> = e', and e kf — e — uj, e k > = + uj, respectively, 
cf. Fig. [j]. Unlike an effective potential mediated by the 
impurity spin [pp| , the t-matrix here is an operator in 
impurity spin space. As usual, the operator T is defined 
by the series 



00 / 1 V 

n=0 V u / 



(9) 



Our interest is in the retarded t-matrix where the 
energy e is determined by the outgoing electrons, i.e., 
e = £fc f + £fe' + iS. Performing perturbation theory up 
to fifth order in the coupling J, we get for the inelastic 
processes 



(10) 



uj 2 2h p 

{l + pJ[g(e)+g(/)+g*(e - cj)+g*(e' + cj)}} 
+l.s. + O(pJ) 6 

where l.s. means less singular terms in 1/lu. To leading 
order in J, this result agrees with the findings of KG. The 
less singular terms omitted arise from energy denomina- 
tors that include an intermediate electronic energy to be 
summed over, leading at most to a logarithmic singular- 
ity for u> — » 0. Further, the auxiliary function g(e) is 
given by 



9(e) 



de' 



/(O 



P(e') 



e - e' + iS 



(11) 



where /5(e) is the normalized density of states of the elec- 
trons with p(0) = 1. 

To go beyond perturbation theory we consider the t- 
matrix in operator form (^). We are interested in the 
most divergent terms for lo — > and therefore search for 
energy denominators of the form l/(e — Hq) = ±1/lu. 
Since e = e kf + = e k - uj + e k , + uj = e k + e k > , 
the intermediate energy must include both e kf and e k i 
or tk' s an d £fe, i.e., two unperturbed electron lines must 
lead to this intermediate state. Therefore, in leading or- 
der in 1/uj we may decompose the two-particle t-matrix 
into two single-particle t-matrices, one acting on the k— 
quasiparticle and the other on the fc'-quasiparticle, cf. 
Fig. 0, 



1 



(12) 



^kajk' a' —*kf a f ^k'^-a'^ T-ha— *h fa j Tk' a' ^-k'^a'^ 

+ ^k'a'^k'a' —Tka^k f a f + 1.S 



Here, the operator character of the t-matrices in impu- 
rity spin space plays an essential role since all commuting 
terms cancel. 
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e'+gd, <f, 



FIG. 1. Decomposition of the two-particle t-matrix into 
two single-particle t-matrices. 

To proceed we follow Zittartz (l7j and decompose the 
single-particle t-matrix into a non-spin flip and a spin 
flip amplitude 



(13) 



Here, Is means the identity in impurity spin space, and 
tkk f and Tkk f are the t-matrices in the non-spin flip and 
spin flip channel, respectively. Since the identity com- 
mutes with other spin matrices, it does not lead to a 
leading order contribution in l/u>. Inserting the two- 
particle t-matrix ( fL?] ) with ( |l3| ) into Eq. (Q) we get 



1 7T 



K(uj,e,e') = — — C imp S(S+l)p 3 \T(e)T(e')\ 2 +l.s. 



(14) 



Comparing this expression with the perturbative result 
(|l0|), we see that the bare coupling J is replaced by the 
renormalized quantity Tkk f — T ( e k f )- This result JT4| ) 
is of central importance for the analysis below. Note 
that the expression (|l4|) displays an explicit 1/uj 2 depen- 
dence which is crucial for the experimentally observed 
scaling behavior || . The kernel for the effective electron- 
electron scattering depends on the off-shell spin flip part 
of the t-matrix which has a different behavior than the 
on-shell non-spin flip part of the t-matrix responsible 
for the temperature dependence of the resistance. In our 
approach the 1/oJ 2 algebraic factor emerges analytically 
from this representation of K(u>, e, e') p2| . 

Zittartz |l7]] has shown that the Nagaoka equations 
lead to a single-particle t-matrix in the spin-flip chan- 
nel of the form 



r(e) 



J 

wr 



(15) 



With Hamann's solution |2l]] of the Nagaoka equations, 
the denominator reads 



(e) = [X{ef + S{S+l){K P jf] 1 ' 2 . 



Here 



X(e) = 1 - S(S + l)(7rpJ) 2 /4 - pJR{e) 
is temperature dependent via the function 



R(e) 



In 



iA 



(16) 



(17) 



(18) 



where A is the electronic bandwidth. R(e) is an equilib- 
rium approximation of the auxiliary function g(e) intro- 
duced in Eq. ( |TT|) . By inserting the solution ( |l5| ) — ( |lg| ) 
into the kernel(ttj) and expanding in J, we can indeed 
recover the perturbative result (pp|). 

Let us now turn to a nonequilibrium situation. Since 
in equilibrium and nonequilibrium the t-matrix expan- 
sions differ only in the occupation probabilities of inter- 
mediate states, we can use the result for the equilibrium 
t-matrix replacing the distribution of intermediate states 
by their (unknown) nonequilibrium form. In the solution 
( |l5| ) — ([IT]) occupation probabilities only effect the func- 
tion i?(e) which in a nonequilibrium system is no longer 
of the form ( |l8| ) but has to be replaced by g(e). Of course, 
then the nonequilibrium single-particle t-matrix in the 
spin flip channel depends via g(e) on the nonequilibrium 
distribution function /(e). The resulting form of the col- 
lision kernel is fully consistent with perturbation theory 
up to fifth order in J and it includes correctly the leading 
logarithmic terms. 

To proceed, the t-matrix and the distribution func- 
tion need to be determined self-consistently. Both quan- 
tities depend on the position x within the wire. We start 
with the initial nonequilibrium distribution function (|^), 
determine the t-matrix from Eqs. ( |l5| ) — (|l7|), and the 
collision kernel from Eq. (|lj). An improved distribution 
function is then obtained from the Boltzmann equation 
(|l|), which gives rise to an improved kernel. This proce- 
dure converges after about 20 iterations. 

First, we compared our results with recent experimen- 
tal data |l3| ] on the energy relaxation in Au wires in pres- 
ence of Fe impurities. In these experiments all parame- 
ters were determined by independent measurements. For 
instance, for sample 2 measured at T = 33mK the diffu- 
sion time to = 1.8ns. The impurity density was deter- 
mined from the temperature dependence of the resistance 



as C\\ 



55 ppm. We set the impurity spin to S = 1/2 



and the Kondo temperature to Tk = IK which is a typ- 
ical value for Fe in Au Further, we used a density 
of states of p = 0.25/siteeV @. 

Pothier et al. || have noticed that the distribution 
functions f(e,V) exhibit scaling behavior so that as a 
function of e/eV all data at a given position x fall on a 
single line. In Fig. || we show our findings for the distri- 
bution function using the parameters of Pierre et al. |l3f| 
for various voltages V = 0.1 — 0.4meV and two positions 
x = 0.25, 0.5. They are compared with the experimental 
data points [ p"3"| and we find excellent agreement without 
adjustable parameters. We emphasize that the result is 
independent of the choice of the bandwidth and is insen- 
sitive to the Kondo temperature as long as T <c Tk ■ For 
other values of x the curves coincide likewise. 

On the other hand, for higher temperatures T > Tk 
where poor man's scaling holds, the kernel varies loga- 
rithmically with the applied voltage V |l2] and the dis- 
tribution function cannot be written as a function of e/eV 
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only. We have also analyzed the data by Pothier et al. 
H on Cu wires and find again good agreement assuming 
a density Ci mp = 14ppm of S = 1/2 impurities with a 
Kondo temperature above 200mK (see Fig. ||) . Small de- 
viations between the theoretical and experimental results 
are likely due to heating in this first experiment [p4[. 
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£/eV e/eV 
FIG. 2. Nonequilibrium distribution function for various 
voltages V and positions x compared with experimental data. 
Due to the scaling behavior of /(e, V), the results fall on a 
single line at each position x. 

To see the effect of the nonequilibrium electronic dis- 
tribution on the renormalization of the exchange inter- 
action, we show in Fig. || for sample 2 of Pierre et al. 
|jj~3f the real part of the single-particle t-matrix in the 
spin flip channel pr(E) for V = 0.2meV at x = 0.5 (left 
panel) and x = 0.25 (right panel). The solid line gives 
the self-consistent solution obtained from the iteration 
explained above, while the dotted line depicts the result 
obtained for the distribution (|J) in absence of energy re- 
laxation. We see that the coupling changes significantly 
with the distribution function. Using only the initial t- 
matrix does not suffice to explain the experimental data. 
Comparing t-matrices for various voltages, we find weak 
dependence between the Fermi points e = ±eV72. This 
gives rise to the scaling behavior of the distribution func- 
tion, while poor man's scaling |L2| implies a significant 
voltage dependence of the scaled data in conflict with 
experiments |13|. 




e/eV e/eV 

FIG. 3. Real part of the t-matrix in the spin-flip channel. 
The solid line shows the self-consistent solution while the dot- 
ted line gives the result using the initial distribution function 
/o. 

To complete our discussion we should have included 
spin dynamics which cuts off the 1/uo 2 divergence of 
K(ui,e,e') for small frequencies |T^| . This is crucial for 
effects that depend strongly on the low frequency limit, 
such as dephasing. However, the distribution function 
/(e) determined from the Boltzmann equation (Q) is al- 
most insensitive to this low frequency cutoff since the 



collision integral remains finite even in absence of a cut- 
off. Hence, for the problem considered here, one can 
disregard the impurity spin relaxation. 

In summary we have determined the effective electron- 
electron collision kernel mediated by magnetic impuri- 
ties for nonequilibrium wires at temperatures below the 
Kondo temperature. We found excellent agreement with 
recent experimental findings |j],[L3|]. In particular, we 
have demonstrated that the distribution function dis- 
plays scaling only in the regime below the Kondo tem- 
perature. 
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